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O ABSTRACT 
O 

^ We present a high spatial resolution Ha survey of 23 cooling flow clusters using 

^ the Maryland Magellan Tunable Filter (MMTF), covering 1-2 orders of magni- 

^ tude in cooling rate, dM/dt, temperature and entropy. We find 8/23 (35%) of 

CN our clusters have complex, filamentary morphologies at Ha, while an additional 

7/23 (30%) have marginally extended or nuclear Ha emission, in general agree- 
ment with previous studies of fine emission in coohng flow cluster BCGs. A weak 
correlation between the integrated near-UV luminosity and the Ha luminosity 
is also found for our complete sample, with a large amount of scatter about the 
expected relation for photoionization by young stars. We detect Ha emission out 
to the X-ray cooling radius, but no further, in several clusters and find a strong 
^ correlation between the Ha luminosity contained in filaments and the X-ray cool- 

ing fiow rate of the cluster, suggesting that the warm ionized gas is linked to the 
^ cooling fiow. Furthermore, we detect a strong enhancement in the cooling prop- 

CN| erties of the ICM coincident with the Ha emission, compared to the surrounding 

^ ICM at the same radius. While the filaments in a few clusters may be entrained 

O by buoyant radio bubbles, in general, the radially-infalling cooling flow model 

OO provides a better explanation for the observed trends. The correlation of the Ha 

Q and X-ray properties suggests that conduction may be important in keeping the 

^ fllaments ionized. The thinness of the fllaments suggests that magnetic flelds are 

an important part of channeling the gas and shielding it from the surrounding 
X hot ICM. 

Subject headings: galaxies: cooling flows - galaxies: clusters - galaxies: elliptical 
and lenticular, cD - galaxies: active - ISM: jets and outflows 
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1. Introduction 

The high densities and low temperatures of the intracluster medium (hereafter ICM) in 
the cores of some galaxy clusters imply massive amounts of gas cooling radiatively out of the 
ICM. Early studies (see review by Fabian 1994) suggested that cooling flows on the order of 
100-1000 M0 yr~^ should be depositing massive amounts of cold gas onto the brightest cluster 
galaxy (hereafter BCG) in the cluster core. The lack of evidence for such large amounts of 
cold molecular gas or young stars is often referred to as the "cooling flow problem" . More 
recent studies suggest that the cooling region is much smaller than previously thought (Peres 
et al. 1998), which implies overall lower estimates on the mass deposition rate. Additionally, 
it is now generally accepted that some form of heating balances the radiative cooling, leading 
to cooling flow rates on the order of 1-10 Mq yr~^. Mechanical and radiative feedback 
from AGN and starbursts (e.g. Peterson & Fabian 2006, Veilleux, Cecil & Bland-Hawthorn 
2005), gas sloshing in the cluster core (e.g. Zuhone et al. 2009), merger shocks (e.g. Randall 
et al. 2002) and conduction from the surrounding ICM (e.g. Voigt & Fabian 2004) have all 
been suggested as possible contributors to ICM heating. However, the relative importance 
of each process remains unclear. 

The presence of optical line-emitting nebulae, found primarily in cooling flow clusters, 
is a clue that cooler (^ 10^ K) gas coexists with the ICM and that cooling flows may not be 
entirely ruled out (Hu et al. 1985, Heckman et al. 1989, Crawford et al. 1999; Jaffe et al. 2005, 
Hatch et al. 2007). However, while several intriguing ideas have been put forward, the origin 
of this gas and the mechanism for heating it are, as yet, unknown. While the presence of 
cool gas could be attributed to the purported cooling flow, one would then naively expect a 
relatively symmetric geometry and a strong connection of optical emission lines with X-ray 
properties. However, this is not generally the case. Crawford et al. (2005) investigate a 
speciflc scenario, the extended Ha fllament in Abell 1795, and suggest as an explanation a 
cooling wake, produced by a symmetric cooling flow falling onto the brightest cluster galaxy 
(hereafter BCG) which has a high relative velocity at the center of the cluster. Oegerle & Hill 
(2001) flnd that ^ 15% of BCGs are moving with a high pecuhar velocity, so an asymmetric 
distribution of warm gas in this fraction of clusters would be relatively straightforward to 
explain in this context. Potentially contributing heating/ionization sources in the filaments 
include (i) the central AGN, (ii) hot young stellar population outside the cD galaxy nucleus, 
(iii) X-rays from the ICM itself, (iv) heat conduction from the ICM to the colder filaments, 
(v) shocks and turbulent mixing layers, and (vi) collisional heating by cosmic rays (Crawford 
et al. 2005; Ferland et al. 2009). 

In star forming regions, the presence of ionized hydrogen is a prime indication of the 
rate and location of star formation. While several cooling flow clusters have been observed 
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to have non-zero star formation rates at their center, based on UV (Hicks & Mushotzky 
2005) and infrared (O'Dea et al. 2008) observations, it is not immediately clear whether 
newly formed stars are fully responsible for the flux and morphology of the observed Ha 
emission. Indeed, there has been very little work done on the spatial correlation between 
the star forming and warm components in cooling flow cluster centers. What has become 
clear from recent studies is that, in general, the measured star formation rate in the BCG 
is consistently less than the spectroscopically determined X-ray cooling rate for the cluster 
by a factor of 3-10 (i.e. Rafferty et al. 2006; O'Dea et al. 2008), implying that not all of the 
coohng gas is turning into stars with a Salpeter (1955) IMF. 

A major problem standing in the way of an explanation of these fllaments is that there 
are still only a small number of known cases with clearly extended emission beyond the 
nucleus of the BCG. This is due to the paucity of narrow-band filters which can be tuned 
to cosmological redshifts. Attempts have been made to quantify the fraction of BCGs which 
exhibit Ha emission using the Sloan Digital Sky Survey (e.g. Edwards et al. 2007). However, 
spectroscopic studies usually only tell us about the very center of the BCG, which can be 
bright in Ha if an AGN is present. What is needed is a deep, wide-field, narrow-band survey 
of clusters of various richness in order to assess the ubiquity, extent, and topology of the 
warm gas. Some progress has been made towards this end by Jaffe et al. (2005), using 
the Taurus Tunable Filter (TTF; Bland-Hawthorn & Jones 1998). The use of a tunable 
filter allows the observer to tune to a large range in wavelengths with a narrow bandwidth 
(^ lOA). The ability to tune the wavelength is particularly useful when observing cooling 
fiow clusters, of which there are very few in the local Universe. 

In the pilot study leading to this survey paper, McDonald & Veilleux (2009; hereafter 
MV09) observed Abell 1795 in Ra with the Maryland Magellan Tunable Filter (hereafter 
MMTF; Veilleux et al. 2010) in Ha. The well-known southeast filament which extends ^ 50 
kpc from the BCG was found to be, in fact, two separate thin strands. These thin filaments 
were resolved into chains of point sources in the far-UV with the Hubble Space Telescope 
(HST) Advanced Camera for Surveys Solar Blind Channel. We suggested that these strands 
of Ha were, in fact, a cooling wake produced by the cooling fiow falling onto a cD galaxy 
which is plowing through the center of the cluster. As the gas cools, it collapses into thinner 
filaments, eventually leading to star formation which will ionize the gas. The thinness of 
the filaments suggests that magnetic fields are likely important to prevent turbulence from 
erasing them, but this is just speculation without any direct evidence of a strong magnetic 
field in the filaments. Although it is likely more complicated than this (we observe higher- 
than- expected [N II] /Ha ratios, and a slight spatial offset between the UV and Ha), this 
scenario can explain the general structure of the Ha complex. 
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In order to have a more complete picture of the extended Ra emission in galaxy clusters, 
we have undertaken a survey of 23 clusters with a wide variety of properties using the 
MMTF. In the following section we describe our sample selection, as well as the acquisition 
and analysis of the data. In §3, we present the results of the analysis and in §4 discuss 
the various implications of these results. Finally, in §5 we present a summary of the results. 
Throughout this paper, we assume the following cosmological values: Hq = 73 km s~^ Mpc~^, 

^matter 0.27, ^vacuum 0.73. 



2. Data Collection and Analysis 

In order to have a more complete picture of cool core clusters, we have undertaken 
a mult i- wavelength survey of 23 cooling flow clusters. From an initial sample size of 232 
clusters, which cover 2 orders of magnitude in dM/dt and with X-ray temperatures ranging 
from 1-12 keV, (full sample from White et al. 1997) we enforced the following cuts: 6 < +35° 
and 0.025 < z < 0.092. These cuts ensured that (1) the cluster was visible with Magellan, 
(2) the appropriate blocking fllters were available to observe redshifted Ha, and (3) the 
cluster angular size did not exceed that of MMTF (^30^). These two cuts removed a total 
of 44 and 62 clusters from the full sample, respectively, leaving a total of 126 observable 
clusters. From this list, we selected 23 clusters to cover the full range of properties, from 
very rich clusters with high cooling rates (e.g., Abell 1795, Abell 2029) to low-density clusters 
with little or no cooling flow (e.g., Ophiuchus, Abell 3376). A summary of these properties 
for our sample clusters is given in Table [TJ As is shown in Table [T} the range in classical 
cooling rates for this sample is from 6.3-431 M© yr~^, which means that, while we cover a 
large range in properties, we are examining only cooling flow clusters. This selection criteria 
should produce a higher fraction of line-emitting BCGs compared to an un-biased cluster 
sample (see Edwards et al. 2007). 

For these 23 clusters, we have data at Ra (MMTF; 23/23), red continuum (MMTF; 
23/23), near-IR (2MASS; 23/23), UV {GALEX, XMM-OM] 21/23) and X-ray {Chandra] 
19/23). For those clusters hosting a radio galaxy, we have also extracted VLA 1.4 GHz 
fluxes from the NVSS (Condon et al. 1998) survey (17/23). The availability of these data 
are summarized in Table [T} 

In the following subsections, we describe in detail the acquisition and reduction of these 
data. 
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Table 1: Sample of 23 cooling flow clusters with MMTF Hct imaging 



Name 


z 


E(B-V) 


Tx 




(1) 


(2) 


(3) 


(4) 


(5) 


Abell 0085 


0.0557 


0.038 


6.5 


108 


Abell 0133 


0.0569 


0.019 


3.5 


110 


Abell 0478 


0.0881 


0.517 


6.8 


736 


Abell 0496 


0.0329 


0.132 


4.8 


134 


Abell 0644 


0.0704 


0.122 


6.5 


136 


Abell 0780 


0.0539 


0.042 


4.7 


222 


Abell 1644 


0.0475 


0.069 


5.1 


12 


Abell 1650 


0.0846 


0.017 


5.1 


122 


Abell 1795 


0.0625 


0.013 


5.3 


321 
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9 

Z.D 
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Abell 2029 


0.0773 


0.040 


7.4 


431 


Abell 2052 


0.0345 


0.037 


3.4 


94 


AKoll 91/19 


n nQfi/i 




1 n 1 
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Abell 2151 


0.0352 


0.043 


2.9 


166 


Abell 3158 


0.0597 


0.015 


5.3 


9.6 


Abell 3376 


0.0597 


0.056 


3.5 


6.3 


Abell 4059 


0.0475 


0.015 






Ophiuchus 


0.0285 


0.588 


8.6 


41 


Sersic 159-03 


0.0580 


0.011 


2.4 


288 


Abell 2580^ 


0.0890 


0.024 


4.3 


95 


Abell 3389^ 


0.0267 


0.076 


2.0 


22 


Abell 0970 


0.0587 


0.055 


4.1 


20 


WBL 360-03"^ 


0.0274 


0.028 


1.8 


10 



(1) - Cluster name 

(2) - NED redshift of BCG (http://nedwww.ipac.caltech.edu) 

(3) - Reddening due to Galactic extinction from Schlegel et al. (1998) 

(4) - Cluster X-ray temperature (keV) from White et al. (1997) 

(5) - Classical cooling rates (Mq yr~^) from White et al. (1997) 
^ - No available Chandra data 

^ - No available Chandra, GALEX or XMM-OM data 
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2.1. Ha: MMTF 

The MMTF has a very narrow bandpass (^5-12A) which can be tuned to any wave- 
length over ^5000-9200A(VeiUeux et al. 2010). Coupled with the exquisite image qual- 
ity at Magellan and the wide field of the Inamori-Magellan Areal Camera & Spectrograph 
(IMACS), this instrument is ideal for detecting emission-line filaments in distant clusters. 
During 2008-09, we observed all 23 clusters at both Ha (A=6563A) and continuum (± 60A), 
for a total of 20 minutes each. If YLa was detected in a 20-minute exposure, it was followed 
up for an additional 40 minutes at each wavelength. The typical image quality for these 
exposures was 0.6 ± 0.2^' 

These data were fully reduced using the MMTF data reduction pipelin^ which performs 
bias subtraction, fiat fielding, sky line removal, cosmic ray removal, astrometric calibration 
and stacking of multiple exposures (following Veilleux et al. 2010; see also Jones, Bland- 
Hawthorn, Shopbell 2002). The continuum image was then PSF and intensity matched to the 
narrow-band images to allow for careful continuum subtraction. The stacked images were cal- 
ibrated using spectrophotometric standards from Oke (1990) and Hamuy et al. (1992, 1994) 
The error associated with our absolute photometric calibrations is ^15%, which is typical for 
tunable filters and spectrographs. Finally, the data were corrected for Galactic extinction, 
following Cardelli et al. (1989) using reddening estimates from Schlegel et al. (1998). We do 
not attempt to correct for intrinsic extinction since the dust content of the optical filaments 
is not well known. All of these procedures are described in detail in Veilleux et al. (2010). 

For systems with complicated morphologies. Ha fiuxes were measured by creating (by 
eye) a region which generously traced the Ha emission and calculating the total signal within 
this region. For more symmetric morphologies, a circular aperture centered on the emission 
peak was used, with the radius chosen to contain all of the obvious emission. Our measured 
Ha fiuxes (see Table A.l) agree well with those in the literature (e.g., Cowie et al. 1983, 
Owen et al. 1995, Edge 2001, Jaffe et al. 2005) for the 12 systems with previously published 
data. In 9 systems we measure slightly higher (^ 1.7x) total fiux, most likely due to our 
better sensitivity and larger field-of-view. In all other cases, the published "Ha" fiux also 
include a contribution from [N II] AA 6548, 6583 and is therefore overestimated. In all cases, 
the Ha emission measured from our data is not contaminated by [N II] emission. 



^ htt p : / / www . astro . umd . edu/ ^veilleux / mmt f / dat ar ed . html 
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2.2. Near-UV: GALEX and XMM-OM 

One possible source of Ha ionization in these clusters is photoionization from hot, young 
stars. In this scenario we would expect to see associated UV emission, perhaps with filamen- 
tary morphology as well. In order to properly assess if this is the case, we have used archival 
GALEX (1350-1750A, 1750-2800A) and XMM-OM (1820-2320A, 2078-2517A, 2595-3215A) 
data for most of the clusters in our sample. The GALEX data for Abell 1795 comes from 
our own program (GALEX GI Cycle 5 Program 31; PI Veilleux). These data were corrected 
for Galactic extinction, following Cardelli et al. (1989) using Galactic reddening estimates 
from Schlegel et al. (1998). 

In order to determine the star formation rate (SFR) from the NUV fiux, we follow the 
prescription described in Kennicutt (1998). This technique assumes that SFR oc with 
no dependence on wavelength (F^ is fiat) over the range of 1500-2800A. In Figure |l| we 
confirm that this assumption is roughly valid for the clusters in our sample and, thus, if dust 
is present it does not follow the galactic reddening law. For 18 clusters, we have GALEX 
near-UV (2267A) data, while the remaining 3 have XMM-OM UVWl (2905 A) data. The 
difference in fiux that we predict between these filters (from Figure [T]) is ^ 1.3 x. Thus, we 
proceed with merging these two sources of UV data with the caveat that UV fiuxes for Abell 
0478, 0496 and 1650 may be -30% too high. 

Accompanying J-band images were also obtained from 2MASS in order to remove the 
expected NUV emission from the old stellar population. Motivated by the techniques de- 
scribed in Hicks et al. (2005, 2010), we attempt to remove the NUV contribution from old 
stellar populations, which we infer from the J-band data. The relationship between the 
GALEX NUV and 2MASS J-band fiuxes follow a powerlaw: 

logio Lnuv = Ci + C2 logio Lj (1) 

Hicks et al. define a non-star forming control sample in order to calibrate this relationship, 
consisting of 17 cluster ellipticals and 22 BCGs in clusters without measurable cooling fiows. 
The BCGs in this control sample cover a similar range in J-band luminosity and redshift to 
the clusters in our sample. Hicks et al. (2010) find Ci = -0.85 ± 0.016 and C2 = 0.89 ± 0.045, 
with a scatter of 0.08, based on fiuxes measured in apertures. Since we are interested 
in detecting extended emission, we instead use a 30^' aperture, which has the net effect of 
changing the zero-point of this relation (due to the NUV- J color gradient). Thus, we adjust 
the zero-point by hand so that it passes through the points with the lowest NUV/ J ratios. 
We find Ci = -0.55 and note that the observed scatter about this fit matches well that quoted 
by Hicks et al. (2010). Fig. |2| shows this fit, along with that found by Hicks et al. (2010), 
for the BCGs in our sample. Using this relation, we scaled the measured J-band fiuxes to 
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obtain the NUV contribution from the old steUar population, and then subtracted this from 
the measured NUV flux. Our value of Ci is slightly larger than that of Hicks et al. (2010), 
due to our use of a larger aperture. We show in the right-hand panel of Fig. [2] that, beyond 
the aperture used by Hicks, the NUV/ J ratio is higher at all radii out to 30'^ Thus, a 
slight increase in Ci is justifled. 

We do not attempt to correct for intrinsic extinction, due to the fact that not all of 
our clusters have reddening estimates, thus the lower envelope in Fig. [2] could also contain 
dusty, star forming galaxies. However, from Crawford et al. (1999), we flnd that the typical 
intrinsic extinction in BCGs with strong Ha emission is E(B-V) = 0.29 ± 0.135. This would 
induce a scatter of 0.5 in the NUV-J relation. Since the measured scatter is much lower 
(0.08), we propose that the majority of the galaxies used in the Hicks et al. (2010) NUV-J 
relation are relatively dust-free. Furthermore, it is unclear how the dust is distributed in the 
BCG (e.g. Lauer et al. 2005) and, thus, any intrinsic extinction correction would require a 
signiflcant amount of guesswork. 



2.3. X-Ray: Chandra 

Archival data from the Chandra X-ray Observatory were retrieved for 19 of our 23 sample 
clusters. These data were reprocessed with CIAO (version 4.1.2) and CALDB (version 4.1.1) 
using the latest time-dependent gain adjustments and maps to create new level 2 event flies. 
Due to the large angular extent of the clusters in our sample, we were required to construct 
blank-sky background event flies, using the ACTS blank-sky background database, to properly 
account for background flux. The new level 2 event flies were cleaned for flares, using the 
Ic-clean routine, by examining the light-curve and removing any spurious bursts in intensity. 
These data cleaning and calibration procedures are all outlined in detail in the CIAO science 
thread^ 

In order to separate any fllaments or interesting morphology from the X-ray halo, we per- 
formed an unsharp masking technique on each image, subtracting a 10'^ Gaussian smoothed 
image from a 1.5^' Gaussian smoothed image. The resulting image highlights any flne struc- 
ture in the X-ray morphology. 

For each cluster, background-subtracted spectra were extracted using dmextract. Up- 
dated response flies were created using mkacisrmf and mkwarf^ following the CIAO science 
threads. Counts were grouped into bins with 20 counts per bin, over the range 0.3 to 11.0 



^ http : / / cxc . harvard . edu / ciao / threads / 
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Fig. 1. — UV SEDs for 15 clusters in our sample with observations at multiple wavelengths. 
All fluxes have been corrected for Galactic extinction. The two GALEX fllters are denoted 
by and G^ while the three XMM-OM fllters are denoted by X\ and The dashed 
lines represent the mean SED slope. These SEDs are sufliciently flat to justify merging 
our 18 GALEX NUV (2267A) and 3 XMM-OM UVWl (2905A) observations into a single 
near-UV sample. 
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log^o Lj [erg s-^] Radius ["] 

Fig. 2.^ Left panel: GALEX NUV (1750 - 2800 A) versus J-band luminosity for the 21 
BCGs in our sample with archival UV data. The errorbars on Lj are smaller than the point 
size in general. The gray dashed line represents the scaling relation determined by Hicks 
et al. (2010) for a 7" aperture, with the scatter in this relation represented by the gray dotted 
line. The black, solid line is our estimate of the contribution to the total NUV luminosity 
from old stellar populations for a 30^' aperture. This line was chosen to fit the several points 
with the lowest NUV/ J ratios, with the assumption that these BCGs have negligible ongoing 
star formation. The dashed lines represent the scatter in this relation determined by Hicks 
et al. (2010). The zero-point difference between this work and Hicks+ is due to the size of 
the aperture considered. By removing this NUV contribution from old stellar populations, 
we can consider only the excess due to young stars. Right panel: Radial intensity profile of 
the BCG in Abell 3158 (a passive, non-star forming BCG) at NUV and J-band. An offset has 
been applied to the NUV intensity to match the J-band luminosity inside of 7" . Beyond this 
radius, denoted by a vertical dashed line, the NUV intensity is consistently higher. Thus, 
by increasing the aperture to 30^^ the NUV/ J ratio should increase. 
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keV. Spectra were extracted in a variety of regions to better understand the relationship 
between the ICM and the Ha emission. These regions are described below, followed by a 
description of the spectral modeling which we apply to these spectra. 

2.3.1. Radial Profiles 

For all clusters, spectra were extracted from a series of concentric annuli. The maximum 
annulus was chosen by eye (typically the largest allowable in the CXO field-of-view) and then 
the interior annuli were chosen iteratively (following Sun et al. 2009) with Tout/^in = 1.25-1.6 
so that each annulus contained a minimum of 10,000 counts. This choice of spacing allows 
for a reliable deprojection and fit to the spectra. The innermost annulus was chosen to have 
an inner radius of 2 pixels, to exclude any contribution from the AGN. The data in this inner 
region of 2-pixel radius is used to identify possible AGN emission. 

2.3.2. Filaments 

In order to properly investigate any possible connection between the ICM and the pres- 
ence of Ha filaments, we also extract X-ray spectra in regions where we detect Ha emission. 
A mask was created from the MMTF images by first registering them to the CXO images 
using the IRAF wregister task. The typical error involved in this process is dependent on 
the relative astrometric calibrations, which was typically a fraction of a pixel, or ^ O.V\ The 
registered image was then smoothed by V using a Gaussian kernel. Finally, all pixels with 
a fiux greater than 3-a above the background and within Rna (defined as the outermost 
radius of Ha emission, measured from the peak of the Ha emission, see §3.1) were included 
in the mask. We extracted spectra in four different regions: (1) < r < Rnai coincident 
with Ha, (2) < r < Rna^ anti-coincident with Ha, (3) O.SRna < R < Rifa, coincident 
with Ha, (4) 0.3Rifc^ < R < R/f^, anti-coincident with Ha. These spectra allow for a direct 
comparison between X-ray derived properties inside and outside of the Ha filaments, both 
near the center and at larger radii. An example of the on-filament extraction regions is 
shown in Figure [3[ for AbeU 1795 and AbeU 0496. 

2. 3. 3. Spectral Modeling 

In order to derive physical quantities from the X-ray spectra, we used the XSPEC 
spectral fitting package (Arnaud 1996). As mentioned above, spectra in radial annuli have 



Fig. 3. — Examples of the on-filament extraction area (shown in red) for AbeU 1795 (left) 
and AbeU 0496 (right). The left-most frame in each pair covers the full radial range, while 
the right-most covers the range O.SRna < r < ^Ha- The off-filament extraction area contains 
all unmarked pixels in the same radial range. 

also been extracted. This allows for the computation of various parameters such as Tx and 
as a function of radius. Since the emission in any given annulus is a combination of emission 
at that radius and also at larger radii seen in projection, we need to deproject the data 
in each radial bin. We do this using two different techniques: direct spectral deproject ion 
(DSDEPROJ; Russel et al. 2008) and the PROJCT routine in XSPEC. We find DSDEPROJ 
to yield more stable radial profiles, in agreement with Russel et al. (2008), who showed 
that PROJCT produces highly fiuctuating radial profiles when applied to a spectrum that 
represents gas at multiple temperatures. Thus the majority of the analysis uses the results 
from this deprojection. However, it should be mentioned that all of the relevant quantities 
which we derive from these profiles are robust to the spectral deprojection method used. 
Along with the spectral deprojection, we used a photoelectric absorption model (PHABS) 
multiplied by a model for thermal bremsstrahlung emission from a hot diffuse gas (MEKAL). 
The elements were assumed to be present in the Solar ratios, as measured by Anders & 
Grevesse (1989). The free parameters for this combination of models are the column density 
of hydrogen, Nh^ the plasma temperature in keV, T^, the metal abundance of the plasma, 
Z, and the normalization of the model, which is directly related to the electron density in 
the plasma, ng. The model normalization {Nmekal) is related to the electron density by: 



where Da is the angular diameter distance to the source, in cm, and dV is the volume of the 
emitting region. We assume that the X-ray halo is spherical, in order to make the volume 
calculation straightforward. In order to determine the electron density inside and outside of 
the ionized filaments, we also need to assume something about their thickness. We discuss 
this problem further in a later section. 

Finally, in order to determine the mass deposition rate in the cooling fiow (M^^ec), we 
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also apply a second set of models to the radially-binned spectra. This new set of models 
consists of the same deprojection and photoelectric absorption mentioned above, multiplied 
by the combination of a thermal plasma (MEKAL) and a second plasma component that 
is cooling (MKCFLOW). The MKCFLOW model has two temperature parameters, high T 
and low T, which we lock to the temperature of the plasma and to 0.1 keV, respectively. 
The choice of the low temperature is arbitrary - we simply need it to be lower than the 
Chandra detection limit 0.5 keV) so that this represents gas that has cooled out of the 
X-ray regime. The metal abundance in the cooling gas is also required to be the same as the 
plasma. Thus, the only free parameter in this new component is the normalization, which 
is precisely the spectrally-determined cooling flow rate, Msp^c- 

2.4. 1.4 GHz Radio: NVSS 

The total 1.4 GHz radio flux at the center of a cluster is a useful diagnostic of AGN 
activity. For the 20 clusters with S > —40° in our sample, we used the 1.4 GHz fluxes 
measured from the NRAO VLA Sky Survey (NVSS; Condon et al. 1998). 

3. Results 

3.1. Warm Ionized Filaments 

This survey is the largest sample of cooling flow clusters with high-resolution, narrow- 
band imaging to date. The success of this survey has been striking, with multiple extended 
filaments seen in 8/23 clusters and slightly extended or nuclear emission seen in an additional 
7/23 clusters. The remaining 8 clusters have no detectable Ha emission. This detection rate 
of ^ 65% is consistent with the value of 7lt.u found by Edwards et al. (2007) for BCGs in 
cooling fiow clusters (for BCGs in both cooling and non-cooling clusters, Edwards et al. find 
the fraction to be ^ 15%). Figure [3] shows the unsharp masked X-ray and archival UV 
data, along with the new MMTF continuum and Ha images, centered on the BCG. Upon 
inspection of these images, the following trends are immediately evident: 

• There is a strong correlation between the X-ray and Ha morphology. Clusters with 
ionized filaments tend to have structure in the X-ray. This structure is brightest at 
lower energies (0.5-2 keV, see Figure [5]), suggesting that the optical filaments live in or 
near cooler X-ray structures. 

• There does not appear to be a correlation between the presence of filaments and the 
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presence of other galaxies near the BCG. In fact, AbeU 2151 and AbeU 2029 each 
have one or more significant companions within the optical radius of the BCG but 
exhibit no detectable Ha emission. Of the 7 clusters with extended Ha filaments, only 
3 (AbeU 0478, AbeU 1644, AbeU 2052) have significant companions within the optical 
radius of the BCG. In a select few cases, the direction of the detected filaments may 
correlate with the position of nearby galaxies (e.g., AbeU 0478, AbeU 1644). 

• The observed warm filaments are never wider than the MMTF PSF (^ 0.6^'), except in 
the case of Sersic 159-03. At a typical redshift of 0.06 (d ^ 250 Mpc), this corresponds 
to an upper limit on the typical filament width of 0.7 kpc. 

• In some cases, where the filaments are very extended (AbeU 1644, AbeU 1795, Ser- 
sic 159-03), there is evidence for matching morphology in the near-UV. This is es- 
pecially true in the case of AbeU 1795, which we have confirmed with much higher 
resolution near-UV images in MV09 and McDonald et al. (in prep; hereafter M+10). 
However, it is challenging to match the NUV and Ha morphologies in general due to 
the mismatch in the PSF widths. 

The properties of the observed Ha filaments are summarized in Table A.l and a detailed 
discussion of individual clusters can be found in the appendix. Based on the properties of 
the Ha emission, the sample can be divided into 3 categories: those with complex, extended 
filaments (type I), those with nuclear emission or only slightly extended (i.e. a few PSF 
widths) emission (type II) and those with no observable emission (type III). We use this 
convention throughout the remainder of our discussion to easily identify trends at other 
wavelengths which could explain both the fiux and morphology of the Ha emission. 

3.2. Star Formation 

A likely explanation for the presence of Ha emission is phot oionizat ion by massive, young 
stars. However, while several studies to date have found star formation in the central cluster 
galaxy, very little work has been done on correlating the observed Ha filaments with young 
stars for a large sample of clusters. As mentioned earlier, archival near-UV data from GALEX 
and XMM-OM exist for 21 of our 23 sample clusters. Figure [6] shows that there is a weak 
correlation (Pearson R = 0.51) between the near-UV and Ha luminosities. This result is not 
surprising since the cooling fiow in these systems is allegedly depositing some cold gas onto 
the BCG and triggering some star formation activity. However, there is considerably more 
scatter than one would expect if the ionization was purely from continuous star formation. 
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Abell 0085 




Fig. 4. — Multi- wavelength data for the 21 clusters in our sample with NUV (GALEX, XMM- 
OM) and Ra (MMTF) data. From left to right the panels are: 1) MMTF red continuum 
image, 2) Unsharp masked CXO X-ray image, 3) MMTF continuum-subtracted Ra image, 
4) GALEX/XMM-OM NUV image. The horizontal scale bar in all panels represents 20 kpc. 
The X-ray and red continuum images are on the same scale, and the Ha and NUV data are 
on the same zoomed-in scale. The square region in the X-ray panels represents the field of 
view for the zoomed-in Ha and NUV panels. 
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Abell 0780 (Hydra A) 




4. — Continued. 
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Abell 1837 




4. — Continued. 
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Abell 2580 




Fig. 4. — Continued. 
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suggesting that modifications to the star formation law (i.e. altering the IMF), or alternative 
physical mechanisms entirely, may be needed to explain the observed NUV/Ha ratios. 

In Figure [6| we also show the NUV/Ha ratios for several different processes, including 
continuous star formation (described by Kennicutt 1998). Using models from Starburst99 
(Leitherer et al. 1999) we age a burst of star formation from 1-20 Myr in order to see the 
evolution of the NUV and Ha fiux. These models assume a metallicity of OAZq and a 
Kroupa IMF (Kroupa 2001) with a = 3.3. The NUV/Ha ratio is quite sensitive to the IMF 
assumed, but the overall trends remain the same. As the massive stars die off, the UV/Ha 
ratio increases rapidly. This scenario is promising, since several clusters have NUV/Ha 
ratios above the value predicted by Kennicutt et al. (1998) and there are very few physical 
mechanisms that can produce such ratios. In this scenario, star formation took place in a 
burst as gas crosses Rcooi and the NUV/Ha ratio is a direct indicator of the time elapsed 
since this burst. If this is indeed the reason for the high UV/Ha ratios, we would expect 
higher resolution UV data to show age gradients along the length of any filaments. This 
trend was not seen in Abell 1795 (McDonald & Veilleux 2009), although there was a slight 
offset between the UV and Ha emission which made quantifying a gradient difficult. An 
alternative explanation for the high NUV/Ha ratios is the dependence of this quantity as 
a function of Ha surface density (Meurer et al. 2009). For star formation in regions with 
low Ha surface brightness (i.e. UGCA44) the UV/Ha ratio is a factor of ^5x higher than 
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Fig. 5. — False color X-ray images of three clusters in our sample. From left to right, these 
are AbeU 0780 (Hydra A), AbeU 1795 and AbeU 2052. The colors are 0.5-2 keV (red), 2-4 
keV (green) and 4-8 keV (blue). The asymmetry in the X-ray ICM is clearly associated with 
the coolest (0.5-2 keV) gas. 

for a region with higher surface brightness (i.e. NGC 1566). Thus, both scenarios of an 
aging burst of star formation and star formation in low surface density environments can 
yield the scatter above the Kennicutt relation that we observe. A third potential source of 
high UV/Ha ratios is intrinsic reddening due to dust in the BCG. We show, in Figure [6| 
the effect of correcting for intrinsic reddening, using extinction measurements from Crawford 
et al. (1999) and Jaffe et al. (2005). Applying this correction tends to move points further 
above the Kennicutt relation, suggesting that the high UV/Ha ratios we measure are not 
due to dust in the BCG. 

Figure |6] presents the global NUV and YLa luminosities, since the poor spatial resolution 
of the GALEX / XMM-OM data does not allow us to discriminate between filament and 
nuclear emission. In MV09, however, we show that the UV/Ha ratio in the nucleus of 
AbeU 1795 is quite low {hpuv /^Ha ^ 7x10"-^^ Hz~-^). In M+10 (in prep), we show that this 
is typical for AGN (central point source) in BCGs based on HST far-UV imaging of cooling 
fiow cluster cores. This means that the effect of including the nuclear contribution in the 
luminosity calculation would move the data below the Kennicutt relation. 

Finally, we show the UV/Ha ratio expected from shock heating (Dopita & Sutherland 
1996). Again, this produces low ratios and thus could only be a strong contributer in very 
few of the clusters. Most notably. Figure [6] shows that the NUV/Ha ratio for AbeU 2052 is 
consistent with shock heating, and we observe that the Ha emission traces the outer shell 
of a cavity carved in the ICM by a radio jet seen at 1.4GHz (see Figure [i]). This is strong 
evidence for the Ha-emitting gas to be shock heated by an AGN outflow in this object. 
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Fig. 6. — Global GALEX NUV versus MMTF Ra luminosities for the 21 clusters in our 
sample with archival GALEX data. The dashed blue line represents ongoing star formation, 
described by Kennicutt (1998). The dotted black lines describe the evolution of the NUV/Ha 
ratio as a starburst ages from 1-20 Myr (Starburst99; Leitherer et al. 1999). The shaded 
grey area shows the range of NUV/Ha ratios expected from shocks (Dopita & Sutherland 
1996), while the shaded red area shows the expected NUV/Ha ratio for the AGN, based on 
HST FUV imaging of cluster cores (M+10 in prep). The orange area represents the overlap 
between these two regimes. The point color refers to the classification described in §3.1: 
type I (blue), type II (green) and type III (red). The dotted lines/circles show the eflFect of 
correcting for intrinsic extinction, with reddening estimates taken from Crawford et al. (1999) 
and Jaffe et al. (2005). The large deviations from ongoing star formation (Kennicutt 1998), 
suggest that non-stellar processes may be responsible for some of the observed Ha emission. 
Contamination from the nucleus likely plays an important role in the observed scatter. 
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3.3. X-Ray Profiles 

Figure [7| shows the radial profiles for various parameters derived from the X-ray data: 
temperature (kT), density (ug), metal abundance, specific entropy (K = kT-n^^^^ keV cm^), 
cooling time {tcooi = 10^-(K/10 keV cm^)^/^-(kT/5 keV)"-^ Gyr) and mass deposition rate 
(dM/dt). There appears to be very little difference in the shape or absolute normalization of 
the density profile between clusters with and without optical filaments. Likewise, while there 
is a large spread in the values, there appears to be very little correlation between the radial 
abundance profile and the presence of optical emission. Interestingly, there appears to be a 
very strong correlation between the absolute level of the temperature profile and our filament 
type classification. That is, for clusters with extended, bright optical filaments (type I), the 
temperature is lower at all radii. For those with nuclear or no emission (types II or III), 
the temperature is higher at all radii. Since the presence of Ha filaments has been shown 
to be enhanced in cooling fiow clusters (Edwards et al 2007), one would expect the core 
temperature to be correlated with the presence of filaments. We find that this correlation 
extends from the core to much larger radius, such that the temperature at all radii from the 
very center out to the cooling radius is lower in clusters with extended Ha filaments. 

The motivation for producing these radial profiles was to look for a correlation between 
the entropy profile and the presence of optical filaments, which was alluded to by Donahue 
et al. (2006). However, while previous authors (e.g. Cavagnolo et al. 2008) have found that 
the central entropy. Kg, is much lower in clusters with Ha emission, we find that the entropy 
is lower everywhere interior to the cooling radius for clusters with extended Ha emission. 
Naively, one would expect strong AGN feedback and cluster-cluster mergers to raise both 
the local temperature and entropy, thus impeding the formation of Ha filaments. However, 
simulations have found that, under the correct circumstances, AGN feedback can lead to the 
formation of cold filaments as radio-blown bubbles rise to large radii (e.g. Fabian et al. 2003, 
Revaz et al. 2008, Pope et al. 2010). 

We observe that both the temperature and entropy are lower in the inner ^ 100 kpc for 
clusters with observable Ha filaments (Figure [7|. In order to quantitatively compare this 
phenomenon to the presence of ionized filaments, we consider quantities derived with this 
scale in mind. Figure [8] shows the comparison of the NUV and Ha luminosities with the 
average temperature inside of 100 kpc (kTioo), the specific entropy at 30 kpc (K30) and the 
mass deposition rate (dM/dt). The UV and Ha luminosities are shown in terms of their 
inferred star formation rates, derived using Kennicutt (1998), so that the comparison to 
dM/dt is straightforward. However, these values are linearly proportional to the luminosity 
and so any observed correlation is independent of this choice of scaling. 

In Figure [8} we see a strong correlation between the UV or Ha luminosity and both 
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10 100 10 100 

Radius [kpc] Radius [kpc] 

Fig. 7. — X-ray derived radial profiles for abundance, electron density (ng), tem- 
perature (kT), specific entropy (K = kT-n^^^^ keV cm^), cooling time {tcooi = 
10^-(K/10 keV cm2)^/2-(kT/5 keV)"^ Gyr) and mass deposition rate (dM/dt). The dashed 
line in the right middle panel represents a cooling time of 5 Gyr. The color scheme here 
is the same as in Figure [6| The lower two panels show clearly that the clusters with YLa 
filaments (blue) have lower temperature and entropy everywhere than those with nuclear or 
no Ha emission (green and red respectively). 
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Fig. 8. — UV and Ha determined star formation rates, based on Kennicutt (1998), against 
the average temperature in the inner 100 kpc (kTioo), specific entropy at 30 kpc (K30), and 
the integrated mass deposition rate (dM/dt). The color scheme here is the same as in Figure 
6| The SFRs in the lower panels have had the nuclear contribution removed (SFR*). The 
open star in the lower panels refers to Perseus A (Conselice et al. 2001, Sanders et al. 2004). 
The four points with dM/dt = 0.3 Mq yr~^ have no measurable cooling flow and are plotted 
as upper limits. The diagonal dashed line represents the limit where all of the cooling X-ray 
gas turns into stars, while the dotted line is the case where all of the cooling X-ray gas is 
made up of hydrogen which recombines only once (Fabian et al. 1984). A strong correlation 
is seen between the Ha luminosity in fllaments, and both the mass deposition rate and the 
average cluster entropy. 
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the entropy and dM/dt over three orders of magnitude. Previous works (e.g. Cavagnolo 
et al. 2008) have pointed out that systems with lower core entropy tend to have UV/Ha 
emission, however we show that the entropy outside of the core is, potentiaUy, an even 
better tracer of Ha emission. We prefer this diagnostic since several of our clusters have 
relatively low surface brightness X-ray emission, making it difficult to constrain the entropy 
in the core where the signal-to-noise is low. The lower entropy systems tend to have brighter 
Ha filaments, with Abell 1795, Sersic 159-03 and Perseus A occupying the lowest-entropy 
and highest Ha luminosity regime. The correlation between dM/dt and the Ha luminosity 
suggests a link between the cooling fiow and either the presence or ionization of the warm 
gas. The observed Ha luminosity is ^1-2 orders of magnitude too high for it to be due to 
recombination from the hot X-ray gas, assuming a single Ha photon per recombination. On 
the other hand, even if all of the Ha is due to star formation, there is not enough ongoing 
star formation to account for all of the gas in the cooling fiow. However, one should be 
cautious when using the X-ray derived dM/dt since they do not represent the true mass 
of gas cooling out of the X-ray, since it is still unclear how cooling proceeds below soft 
X-ray energies (Peterson & Fabian 2006). Note, nevertheless, that there appears to be a 
temperature threshold of ^ 5 keV for ionized filaments, below which nearly all clusters in 
our sample have extended Ha emisison. It is unclear at present whether there is a lower 
temperature threshold below ^ 2 keV, which is the domain of galaxy groups. 

From the radial X-ray profiles and the 2-dimensional Ha data, we are also able to derive 
two interesting scale radii. These are: Rcooh the radius at which the cooling time of the X-ray 
gas is equal to 5 Gyr, and R/f^, the maximum extent of the Ha filaments. The relationship 
between these radii is shown in Figure [9} While there is a large observed range in the radius 
of the Ha emission, there appears to be a hard edge at Rcooi ~ we never observe Ha emission 
beyond the radius where the X-ray cooling time is > 5 Gyr. This is a strong indication 
that the Ha emission is intimately linked to the X-ray cooling fiow. The fact that the data 
appear to prefer a certain timescale (5 Gyr) for the cooling radius is a new and puzzling 
discovery. Interestingly, recent studies have shown that beyond z ^ 0.5 the frequency of 
clusters with cool cores is significantly lower than in the local Universe (see e.g. Vikhlinin 
et al. 2007, Santos et al. 2008). This redshift corresponds to a lookback time of 5 Gyr, 
which provides further motivation for our choice of cooling radius. Thus, if cooling fiows 
were unable to begin until 5 Gyr ago as the evidence suggests, we would expect the radius 
at which tcooi = 5 Gyr to be a natural boundary for the Ha filaments. 

Additionally, if we adopt the idea of stars forming in bursts, as indicated by the location 
of the points in Figure [6| our data show that the emission from the youngest stars tends to 
be much more extended, while the emission from stars formed in a burst ^ 10 Myr ago tend 
to be concentrated in a nuclear region, lending some support to a radially-infalling scenario. 
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Fig. 9. — Correlation and distributions of Rcoo/, the radius at which the coohng time of the 
ICM reaches 5Gyr; and Rij^, the largest radius at which we detect Ha emission. For the 
cases with no detected Ha emission an upper limit of the seeing FWHM has been used. The 
dashed line refers to the one-to-one case, while the point colors are consistent with Figure [6| 
The shapes of the symbols refer to the NUV/Ha ratio: circles have NUV/Ha < 10~-^^-^ Hz~-^, 
while triangles have NUV/Ha >10~^^"^ Hz~^. Note that, in the burst formation scenario 
this corresponds to age, with circles representing younger and triangles older stars. There 
appears to be an upper limit on the radius of Ha filaments, corresponding to the cooling 
radius. The black point represents a lower limit estimate on the cooling radius of Perseus A 
(Fabian et al. 2000, Conselice et al. 2001). The effect of decreasing the cooling time used to 
defined Rcoo/ is to move points down in this plot. 



-27- 



Csl 
I 

U 
CD 
if) 
U 



u 
T 

CO 
CD 



CN 
I 

u 

CD 
CO 
U 



CD 



I 

10-1^ 

10-18 
10-19 



T — I — r 



n — I — r 



T — r 



Abell 1795 == 



J \ L 



1^ 



J \ L 



L_ L 



^ I ' [ ^ 
Sersic 159r03 ^ 




1 x^^f 







20 40 
r ["] 







10 



20 
r ["] 



30 



10-15 

D 10-16 

M 

I 

E 10-17 

10-18 
10-19 



Abell 2{)52 =| 



:X « 



,>*H< 



• I 



X 

X X 



X 

J >L. 



Abell 0496 ^ 



X X 



_LJ L 







10 20 







20 
r ["] 



40 



Fig. 10. — Ha surface brightness profiles for the three clusters with RRa/^cooi ^ 1 and 
one cluster with Rna « ^cooi (Abell 0496). The black crosses represent the azimuthally 
averaged profile, while the blue circles are the surface brightness along a ^ 5^' wide cut, 
oriented along the most extended filament. The vertical dashed line represents the location 
of Rifa, as measured by eye. The sharp drop-off of surface brightness at Rna suggests that 
this is truly the edge of the filament, independent of the depth of our exposure. 
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Since we are claiming a link between the cooling radius and the maximum radius of Ha 
emission, it is important to assess whether we are truly seeing all of the Ha emission, or if a 
deeper exposure would yield RHa/^cooi > 1. In Figure 10, the YLa surface brightness profiles 
of the three clusters with Rna/^cooi ^ 1 are shown. At Rna^ there is a clear drop-off in the 
surface brightness in all three systems, suggesting that this is a real limit, not a function 
of our exposure depth. For clusters with Rna < 0.5Rcooi (e.g. Abell 0496, Figure 10, the 
edge of the filament blends into the background much more smoothly, suggesting that these 
filaments may in fact extend even further. 



3.4. Properties of X-Ray Gas On and Off of Filament 

Using smoothed Ha images as masks, we were able to extract CXO spectra on and sur- 
rounding the warm ionized filaments, as described in §2.3.2. The off-filament spectra were 
fit with a single-temperature model, while two cases were considered for the on-filament 
spectra. In order to quantify the so-called "iron bias" (Bohringer et al. 2009), which re- 
sults from fitting a multi-temperature spectrum with a single-temperature model, we try 
fitting both one-temperature and two-temperature models to the on-filament spectra. A 
two-temperature model is clearly motivated if we assume that the ICM conditions are dif- 
ferent on-filament, since it allows us to account for the ICM seen in projection in front of 
and behind the filaments. The S/N of the on-filament spectra are typically insufficient to 
constrain the temperatures, abundances and relative normalization in a two-temperature 
model, so we opt to freeze the temperature and abundance for one of the two models to the 
off-filament conditions. 

As we mentioned in §2.3, measuring the electron density on-filament requires an assump- 
tion about the geometry of the extraction region. For circular annuli, this is straightforward 
since we assume a spherical cluster. Since we can not be certain about the geometry of the 
filaments, we consider two limiting cases: 1) the filaments are thin tubes and extend no fur- 
ther than the width of our MMTF PSF 1 kpc) into the plane of the sky; 2) the filaments 
are sheets seen edge-on with length equal to the observed size of the cluster. For the first 
case we model the spectrum with a two-temperature plasma to represent the thin filament 
embedded in the surrounding ICM, while the second case is modeled as a single-temperature 
slab of gas. The fact that we only see filaments that are unresolved in width suggests that 
the filaments are indeed thin and not sheets seen in projection. However we proceed un- 
der the assumption that both geometries are reasonable, thus bounding the "true" solution 
somewhere in between. 

Fortunately, the assumption of geometry does not affect the calculation of the abun- 
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dance or temperature. The on-versus-off filament temperature and abundance ratios are 



sliown in Figure [11} for both single-temperature and two-temperature models. Immediately 
obvious from this figure is the fact that the filaments tend to be cooler than the immediate 
surroundings - the typical temperature in the filament is ^ 50% of that in the nearby gas. 
The fact that the X-ray gas near the Ha filaments is significantly cooler than the surrounding 
ICM is further evidence that the Ha filaments may be taking part in the cooling fiow. Addi- 
tionally, the on-filament metallicity tends to be significantly lower, reaching values less than 
^ 40% of the off-filament values for about half of our systems. This does not appear to be 
the case for the 2-temperature model applied to the "outer filament" spectrum, suggesting 
that the low on-filament metallicity may be a result of the iron bias (Bohringer et al. 2009). 
However, this is also the lowest S/N of the four scenarios considered. Thus, deeper data are 
required in order to say for certain what the true abundance gradient between the on- and 
off-filament spectra is. There are very few physical processes that decrease metallicity so, if 
the on-filament abundance is in fact lower with respect to the surrounding ICM, it is more 
likely that the low-metallicity gas has been transported from elsewhere. 

A hint as to the origin of this low-metallicity gas can be found if we look at the 5 
clusters with the largest on-off filament metal abundance contrast. With the exception of 
Abell 2052, all of these clusters have very strong radial abundance gradients. On the other 
hand, those clusters which do not have a strong on-off filament metal abundance contrast 
have relatively fiat radial abundance gradients. To further drive this home, we show in 
Figure [12] the relationship between the on- and off-filament abundance and the abundance 
in four regimes: (1) r ^ 0, (2) r < icooh (3) r ^ icooi and, (4) r > icooi- This figure shows that, 
in general, the metallicity in the filaments matches best with the metallicity at TcooI- For 
those four clusters with strong on-off filament abundance gradients (excluding Abell 2052), 
the on-filament abundance matches almost exactly the local abundance at TcooI- For the 
remaining four clusters, the on-filament abundance looks roughly the same as the abundance 
at all radii, due to the relatively fiat abundance gradient. If we consider a two-temperature 
model for the on-filament spectrum, the correlation between the on-filament metallicity and 
the metallicity at the cooling radius becomes considerably less significant. Again, more 
data is needed in order to achieve a high enough on-filament S/N in order to constrain the 
abundance in a two-temperature model. A direct link between the on-filament metallicity 
and the metallicity outside of the cooling radius would further strengthen the connection 
between the Ha filaments and the cooling fiow. Furthermore, it would offer some evidence 
against the common hypothesis that the gas in the Ha filaments comes from the central 
region, either from radio-blown bubbles (i.e. Revaz et al. 2008) or from gas sloshing (i.e. 
Johnson et al. 2010). While in some cases the match between the core and and on-filament 
abundances is good, it is overall worse than the match to the abundance at Rcooi- However, 
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Fig. 11. — Upper left panel: Ratio of metallicity and temperatures of the X-ray gas coincident 
with Ra filaments to that of the surrounding ICM for a single-temperature model. The 
dashed lines represent equal temperature and abundances inside and outside of the filaments. 
The filaments tend to have systematically lower temperatures and abundances in all cases 
with the exception of Abell 0780. Upper right panel: Similar to the left panel, except that 
only the outer 70% of the filament in radius is considered (see §2.3.2). Lower panels: Similar 
to the upper panels, except using a 2-temperature model to fit the on-filament spectrum. 
Note that this results in significantly lower on-filament temperatures and overall larger errors 
in both fitting parameters. 
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Fig. 12. — Upper panels: Metallicity of X-ray gas coincident with the Ra filaments as a 
function of the metallicity in various radial bins. The number in each panel is the standard 
deviation around the line of equality. The metallicity of the in-filament gas matches more 
closely the metallicity of the ICM at the cooling radius than in any other radial region. The 
lowest point in the upper panels, which is an outlier in all four panels, is Abell 2052. For 
a discussion of why this cluster is an outlier, see the section on indixidual clusters in the 
Appendix. Middle panels: Same as upper panels, but now considering a two-temperature 
fit to the on-filament spectrum. The additional free parameters yields a poorly constrained 
in-filament abundance. Again, the abundance at the cooling radius appears to match best 
the in-filament abundance, but the difference is now only marginal. Lower panels: Same as 
above, but for the X-ray gas outside of the Ra filaments. 
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Fig. 13. — Upper panels: Distribution of temperature, electron density, entropy and cooling 
time ratios in and out of the Ha filaments. The two line types bracket the extreme cases of 
the filament geometry: solid black lines are the thin filament case, which is modeled with a 
two-temperature plasma, while the red dashed line are the case of single-temperature sheets 
seen edge-on. The in-filament gas shares similar properties with the cooling fiow, namely 
that the X-ray gas has a cooling time of ^20% that of the surrounding ICM. Lower panels: 
Similar to the above panels, but now considering only the outer 70% of the filaments in 
radius. 
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there are two biases which are hkely conspiring to increase the uncertainty in the core 
abundance measurement. First, we expect the core abundance to be most affected by the 
iron bias, implying that the core abundances which we quote are hkely lower limits. At the 
same time, our coarse binning in the central region (chose to ensure a high S/N deprojection) 
means that any strong dip in the central abundance will be smeared out, leading to a slightly 
higher measured abundance. As an example, for Abell 4059 we find a dip in the abundance 
profile in the center, with a central value of ^ 1.0, while Reynolds et al. (2008) quote a value 
of ^ 0.3 due to their higher resolution. Using a 2-component model, Reynolds et al. find 
that the core abundance dip changes to an excess, with a central value of ^ 2.5. Thus, it is 
clear that the value of the central abundance is highly dependent on the model and binning 
chosen, and is therefore uncertain. 

In addition to the X-ray temperature and metallicity in the filaments, we can also 
measure various additional properties such as electron density, entropy and cooling time, as 



shown in Figure 13, The distribution of these properties offers further evidence for a link 
between the X-ray cooling fiow and the Ha emission. Assuming the thin-filament geometry, 
the density inside of the filament is typically a factor of a few higher than the surrounding 
ICM. Coupled with the temperature decrease mentioned above, this yields drastically lowered 
entropy and cooling times inside of the filament. With cooling times ranging from 10%-40% 
of the surrounding ICM, these regions are, by definition, part of the cooling fiow. It should 
also be noted that, even if we assume that the geometry is that of single-temperature sheets 
seen in projection, the cooling time and the entropy are lower overall inside of the filaments. 

The overall results from extracting X-ray spectra coincident and anti-coincident with 
the observed Ha filaments tell an interesting story. The Ha filaments appear to reside in a 
portion of the ICM with higher densities that is cooling more rapidly than its surroundings, 
resulting in lower temperatures. The metallicity inside of the filaments resembles more 
closely the metallicity at the cooling radius than the adjacent off-filament gas, albeit with 
a considerable amount of uncertainty. These two results, coupled with the fact that the 
observed Ha emission seems to know about the cooling radius {Rna ^ Rcooi)^ suggests that 
the observed Ha filaments are coupled to the cooling fiow. 



3.5. Properties of BCGs 

As a method of probing whether or not these clusters harbor an AGN in their central 
galaxies, we consider the 1.4 GHz radio power as well as the presence of an X-ray point 
source in the cluster core. The distribution of radio luminosities for our sample is shown in 



Figure 14, We note that all clusters with filaments have a significant radio detection, while 
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Fig. 14. — 1.4 GHz radio luminosity for all clusters detected in the NVSS, as a function of 
the filament type described in §3.1. The red star-shaped points represent the two clusters 
with an X-ray point source detected in the center, while the filled black circles represent 
those clusters with no detected X-ray point source. There appears to be little correlation 
with the radio power and the presence of optical filaments, as 2 of the 6 brightest clusters 
in radio have no detectable Ha emission whatsoever. 
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Fig. 15. — Ha flux contained in fllaments (nucleus removed) versus the total BCG K'-band 
flux from 2MASS. There is no obvious correlation between the brightness of the BCG and 
the amount of warm ionized hydrogen in fllaments. 



only 5/15 of the clusters without fllaments have a signiflcant radio detection. Thus, the 
presence of radio emission is not predictive of whether there will be extended Ha fllaments. 
The two clusters with detected X-ray point sources in the BCG, Abell 0780 (Hydra A) and 
Abell 2052 have large cavities in their X-ray halos, suggesting that the AGN is strongly 
influencing the surrounding medium. The Ha in these two clusters is very likely linked to 
this activity. If we consider the remaining clusters, the correlation between radio power and 
the presence of Ha fllaments becomes much weaker, suggesting that the two phenomena are 
not directly tied to one another. Instead, whatever weak correlation that is seen could be 
due to the fact that cool core clusters, which tend to have Ha fllaments (Figure [s]), will also 
have a more ready supply of fuel for the central AGN, leading to increased radio feedback. 
Further discussion of the implication of these, and other, results follow in the next section. 

We can also consider the relationship between the BCG itself and the presence of optical 
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filaments. Figure [15] sliows tlie total K'-band flux for the BCG versus the total Ha flux 
contained in fllaments. The K'-band flux is a good approximation to the total stellar mass, 
with only a very small deviation with galaxy color (Bell & de Jong 2001). There is no 
obvious correlation between the the galaxy brightness (or stellar mass, approximately) and 
the presence of Ha fllaments, suggesting that the observed fllaments are related instead to 
the X-ray properties of the cluster core. 



4. Discussion: Origin of the Ha Filaments 

In the previous sections we have provided several new clues to the origin of the observed 
Ha fllaments: 

• In clusters with fllamentary Ha emission, we tend to see structure in the cooler (0.5- 
2.0 keV) X-ray gas. If there is no structure in the X-ray gas, there is typically no 
accompanying Ha emission (Figures [31 [5|). 



• There is a strong correlation between the Ha flux in fllaments and both the X-ray 
determined cooling flow rate, dM/dt, and cluster entropy over 3 orders of magnitude 
(Figure [8]). 



• The extent of the Ha fllaments never exceeds the cooling radius. This appears to be a 
hard limit to the radius of these fllaments (Figure [9]). 

• The metallicity of the X-ray gas coincident with the Ha fllaments is often lower than the 
non-coincident ICM at the same radii. The on-fllament X-ray metallicity is consistent 



with that measured near the cooling radius or beyond (Figure 12) 



• The temperature and cooling time of the X-ray gas coincident with the Ha fllaments 



is signiflcantly lower than that of the surrounding ICM (Figure 13). 



These results suggest a direct link between the hot ICM and the warm ionized fllaments. 
However, before we discuss their implications, we would be remiss if we did not mention the 
spectacular Ha fllaments detected in the nearby cooling flow cluster, Perseus. A great deal 
of literature has been devoted to the study of these fllaments, which are seen in molecular 
gas (Salome et al. 2008), warm ionized hydrogen (Conselice et al. 2001) and X-ray (Fabian 
et al. 2006). However, while the observed fllaments are indeed impressive, we would argue 
that they are "typical" among what we classify as our type I clusters. In Figure [16] we 
show the effect of placing Perseus A at larger redshift, and simulating a 20-minute exposure 
using the MMTF under typical observing conditions. By a redshift of 0.06, which is roughly 
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the mean redshift of our sample, the majority of the fine structure has been washed out, 
leaving only a small number of detectable filaments. These redshifted images look similar 
in morphology to our detections of filaments in the 7 clusters which we label as type 1. 
Additionally, we show in Figure|8]that Perseus A occupies a similar space to the other clusters 
with filamentary Ha in temperature, cooling rate and temperature. The total luminosity in 
Ha is significantly higher than in all other clusters, however this is probably again partially 
due to its proximity: if a filament completely disappears into the background we will choose 
a smaller aperture and, thus, measure a smaller total luminosity. 



z=0.0176 224" 


z=0.03 131" 


z=0.06 68" 


z=0.1 39" 



















Fig. 16. — Image of Perseus A from Conselice et al. (2001). The image quality has been 
degraded to simulate the effect of observing this cluster at varying distances. The redshift of 
observation is in the upper left corner, while the apparent size of the image in the horizontal 
direction is in the upper right. By the mean redshift of our sample, z ^ 0.06, most of the 
spectacular filaments have been lost, leaving only a few unresolved structures. 

We finish by discussing the results summarized above in the context of several current 
theories which attempt to explain the presence of the observed Ha filaments: buoyant radio 
bubbles, gas sloshing in the central potential, accretion of gas-rich galaxies, filamentary 
cooling fiows, ICM conduction and magnetic fields. 



4.1. Origin of the Cool Gas 

4.1.1. Buoyant Radio Bubbles 

In an effort to explain how an AGN can deposit energy into the ICM and quench cooling, 
much work has been focused on the evolution of bubbles blown in the ICM by radio jets 
(e.g. Reynolds et al. 2005, Vernaleoet al. 2007, Revaz et al. 2008). Due to the local density 
contrast, these bubbles are buoyant and will rise to larger radius, transporting energy from 



-38- 



small radii to large. While rising, the shape of the bubble evolves and can leave behind 
a trail of cooler gas that morphologically matches the observed Ha filaments in Perseus A 
(Reynolds et al. 2005, Revaz et al. 2008). Radio bubbles have been observed in several 
clusters to date, including some in our sample (e.g. Abell 0780, Abell 2052), which lends 
further support to this scenario. 

This scenario is consistent with the low temperatures seen in the X-ray ICM coincident 
with the Ha filaments and with the fact that all of the clusters with extended filaments 
have non-zero 1.4 GHz fiux. Some clusters (e.g. Abell 0780, Abell 2052) show Ha emission 
along the edges of known radio-blown bubbles lending further support to this scenario. It 
is unclear exactly how the strong correlation between the Ha radius and the cooling radius 
fits in this scenario, although it could be that this relation is tied to the ionization and not 
to the source of the gas. The low ICM abundance coincident with the Ha filaments appears 
to match well with the abundance beyond the cooling radius, but it may also match the 
core abundance, which previous works (e.g., Sanders et al. 2002, 2007) have found to be 
metal-poor in several nearby clusters. 



4' 1-2. Gas Sloshing 

The gravitational disruption of a cluster core by mergers or interactions with other 
clusters or galaxies can lead to the ICM "sloshing" about the central potential, creating 
cold fronts that have been observed in several clusters to date (e.g. Zuhone et al. 2009, 
Johnson et al. 2010). While this is certainly an important mechanism for heat transport 
and is very likely ongoing in several clusters with observed cooling fronts, it does not seem 
to be a plausible general explanation for the presence of warm filaments. The morphology 
of the Ha filaments in our sample tends to be radial, with very few exceptions. The gas 
sloshing scenario predicts cool, low-entropy gas along fronts which are aligned perpendicular 
to the radial direction. The process of sloshing in the cluster core should also smooth out 
any abundance gradients, however we detect possible abundance contrasts between the X- 
ray gas on-filament and off-filament. Finally, the magnitude of the gas sloshing should be 
dependent on the strength of the merger or disruption event. However, we see a clear limit 
to the length of the filaments (Rna < ^cooi)^ which is inconsistent with the sloshing picture. 
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4.1.3. Accretion of Younger Galaxies/ Groups 

While the direction of some of the observed filaments seemingly coincides with the 
location of smaller galaxies (e.g. Abell 0478, Abell 0496, Abell 0780), it is unlikely that the 
warm gas has been stripped from these galaxies. If this were the case, one would expect 
the filaments to extend all the way to the satellite galaxy, in some cases, which is not seen. 
Likewise, we would not expect a single galaxy, or small group of galaxies, to have a strong 
effect on the distribution of the hot X-ray gas. However, much of the structure seen in Ha 
is seen also in the Chandra X-ray images. Additionally, the strong correlation between the 
mass deposition rate and cluster entropy would be unexplained by such a formation method. 

A back-of-the-envelope calculation seems to rule out this mechanism as the source of 
the observed Ha filaments. The typical Ha luminosity that we observe in filaments is ^ 
5x10^^ erg s~^, which corresponds to ^ 3.7x10^^ recombinations per second (assuming case 
B recombination). If we assume that the diffuse ionized medium in the stripped galaxy is 
pure hydrogen and is being stripped from a galaxy like our Milky Way, and then recombining 
only once as it cools in the absence of ionizing radiation, this corresponds to ^ 950 yr"""^ 
of recombining gas. Assuming a total diffuse HII mass of 10^ Mq for a typical massive spiral 
(Lequeux 2005), this implies that all of the stripped gas would recombine in ^ 1 Myr. Thus, 
in order to consistently produce warm ionized filaments via gas stripping, the BCG would 
have to strip the gas disk from a spiral galaxy every few Myr, which is unreasonably short, 
given the typical cluster crossing times of ^ 1 Gyr. Alternatively, if we invoke an additional 
ionization source once the gas has been stripped from the galaxy (e.g. X-ray conduction 
from the ICM), a smaller amount of stripped gas could lead to the observed fluxes of Ha. 
We will revisit this issue in much greater detail in subsequent papers. 

4.1.4' Filamentary Cooling Flows 

In this scenario, the observed Ha fllaments directly trace the cooling flow (i.e. Cowie 
et al. 1980; Fabian et al. 1984; Heckman et al. 1989). The X-ray ICM begins to cool 
dramatically at the cooling radius, while falling inwards at roughly the free-fall speed. As 
the gas cools, small density contrasts in the ICM grow and thin, high-density fllaments 
form. Once sufflciently cool, the gas in these fllaments can be re-ionized via star formation, 
coUisional heating (e.g. Ferland et al. 2009) or drag heating (e.g. Pope et al. 2008). This 
scenario is consistent with Ha fllaments extending to Rcooh with the observed correlation 
between X-ray and Ha properties and with the low metallicity and high cooling rates found 
in the X-ray gas coincident with Ha fllaments. Additionally, the fact that we only see Ha 
fllaments out to a radius where gas can cool in 5 Gyr agrees well with observations that the 
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frequency of cool cores beyond a lookback time of 5 Gyr is significantly lower than the local 
Universe. Thus, the cooling flow scenario is our preferred formation scenario for the optical 
fllaments. 

The thinness of the observed Ha fllaments relative to their extent (Table A.l) is remi- 
niscent of the thin, fllamentary structure seen in hydrodynamic simulations of cool gas flows 
(e.g. Hattori et al. 1995, Pope et al. 2008, Ceverino et al. 2010). Figure 3 from Ceverino 
et al. (2010) shows long (> 50 kpc), thin (< 10 kpc) filaments of cold gas which resemble 
our observations of Ha filaments. While the physical scales are slightly different, they de- 
scribe the process which cool (< 5x10^ K) filaments embedded in a spherical halo of hot 
(> 3x10^ K) gas may form. The formation of these cool filaments is used by the authors 
to explain how star formation proceeds in spiral galaxies at high redshift, but the physical 
processes should be similar for gas coohng from the hot phase to the cool phase, regardless of 
whether we are considering gas streaming onto an isolated spiral galaxy or a central cluster 
galaxy. 

In MV09 we argued that the most plausible scenario for the creation of the long, thin 
filament in Abell 1795 was a cooling wake caused by a cooling flow occurring around the 
moving cD galaxy. This claim was made based on the observation by Oegerle & Hill (2001) 
that Abell 1795 has a signiflcant peculiar velocity relative to the cluster. However, we find 
that the second most asymmetric case, Abell 0496, which also has filaments extending only in 
one direction, has a negligible peculiar velocity. Likewise, several BCGs with no Ha emission 
(e.g., Abell 2029, Abell 2151) or nuclear-only Ha emission (e.g. Abell 2580) have a peculiar 
velocity as high, or higher, than Abell 1795. Thus, while this may play a role in determining 
the direction of the filaments, it is not likely the most important factor in their formation. 

This scenario makes the prediction that the filament material should be in near free- 
fall. The jury is still out on this issue, as several works (i.e. Hatch et al. 2007, Oonk 
et al. 2010) argue against radial infall in Ha filaments, while others (i.e. Lim et al. 2008, 
Wilman et al. 2009) argue for it. Thus, a homogenous sample of kinematic measurements 
is needed to complement this census of the warm and hot medium and provide a strong 
argument either for or against the radially-infalling cooling flow. 

4.2. Evidence for Conduction 

Conduction from the X-ray ICM is a promising scenario which could provide the energy 
necessary to ionize the cooling gas. Nipoti & Binney (2004) predict that, if conduction is 
responsible for ionizing the Ha fllaments, there should be an excess of soft X-rays surrounding 
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the filaments. This is consistent with our findings that the soft X-ray morphology matches 
well with the Ha morphology and that the temperature of the X-ray gas coincident with the 
warm filaments is significantly cooler than the surrounding ICM. In an attempt to address 
the plausibility of this scenario for ionizing the observed Ha filaments in our sample, we 
turn to a study of thermal conduction in 16 cooling-fiow clusters by Voigt & Fabian (2004). 
Of these 16 clusters, there are 5 which overlap with our sample: Abell 0478, Abell 0780, 
Abell 1795, Abell 2029 and Sersic 159-03. Of the five clusters in this overlapping sample, 
four show very extended, filamentary Ha emission. The conductivity inside of the cooling 
radius for these four clusters ranges from Keff/i^s — 0.3-2.3, while for the lone overlapping 
cluster with no Ha emission, the effective conductivity is i^eff/i^s — 0.2 (Voigt & Fabian, 
2004). Estimates of the minimal value required for conduction to to transport heat from 
the hot ICM to cooler regions varies from 0.01 (Churazov et al. 2001) to 0.3 (Narayan & 
Medvedev 2001). 

The fact that we observe cooler X-ray gas coincident with the Ha filaments, coupled 
with the relatively large conductivity measurements in clusters with filaments (albeit with 
a small sample size) suggests that conductivity may be responsible for much of the energy 
required to ionize these filaments. However, future work is needed to expand on the number 
of clusters with effective conduction measurements and to determine if the energetics of the 
warm filaments match what we would expect from ICM conduction. 

4.3. Role of Magnetic Fields 

Given the typically large axial ratio of the Ha filaments (length/thickness ^ 30, accord- 
ing to Table A.l), magnetic fields likely play a role in the formation of these filaments, in one 
or more of several ways: (i) by channeling the infiow of clumps of gas along field lines, giving 
them coherent structure (Fabian et al. 2003), (ii) preventing the hot, turbulent ICM from 
shredding the filaments (Hatch et al. 2007), (iii) to help the growth of thermal instabilities, 
leading to thin high-density filaments (Hattori et al. 1995), and (iv) suppressing conduction 
of heat from the ICM to the coohng gas (Voigt & Fabian 2004). 

High-resolution X-ray spectra of these clusters could constrain the level of turbulence 
in the ICM (e.g., Sanders et al. 2010). These measurements, combined with the observed 
thinness of the Ha filaments and MHD models could constrain the strength of the intracluster 
magnetic field required to sustain such filaments, in a similar manner to Loewenstein & 
Fabian (1990). This would be a nice consistency check on the magnetic field strengths 
derived from rotation measures. 



5. Summary and Future Prospects 



This study represents the highest resolution survey of warm ionized filaments in cooling 
fiow clusters beyond 100 Mpc to date. Complex, filamentary morphologies are seen in 
8/23 (35%) of our clusters, with slightly extended or nuclear emission seen in additional 
7/23 (30%). We find a weak correlation between the total Ra luminosity and the near-UV 
luminosity from young stars, which suggests that photoionization from young stars may play 
a role in producing the observed Ha. UV imaging of higher spatial resolution is needed 
to spatially correlate the Ha and UV emission, and to remove the contribution to the UV 
luminosity from the AGN and/or central starburst. While the Ha and X-ray data support 
the buoyant radio bubble hypothesis as a source of cool gas in some cases, we believe that the 
cooling fiow model does a better job of explaining the properties of the ensemble of clusters. 
We observe correlations between the soft X-ray morphology and the Ha morphology, the 
X-ray cooling rate and the Ha luminosity contained in filaments. We find evidence for a 
correlation between the X-ray metal abundance of gas coincident with the Ha filaments and 
at the cooling radius, however more data are required in order to obtain reliable on-filament 
abundance measurements. The maximum radial extent of the Ha emission and the cooling 
radius appear to be linked, such that Rna < ^cooh where Rcooi is the radius at which the 
X-ray gas can cool in less than 5 Gyr. This radius appears to be physically motivated, as the 
choice of timescale agrees well with the lookback time at which cooling fiow clusters begin 
to emerge. Additionally, we find that the X-ray gas coincident with observed Ha filaments 
is cooling on a much shorter (^ 20%) timescale than the surrounding ICM. These results 
suggest that the Ha filaments trace the X-ray cooling fiow, and that this fiow is initiated at 
the cooling radius. There does not appear to be a strong correlation between the presence of 
filaments and the structural properties of the BCG, suggesting that the Ha filaments are not 
related to any feedback processes from the central galaxy. Based on our results, we propose 
that the gas is cooling out of the X-ray ICM in clumps at roughly the cooling radius and 
collapsing into thin streams as it falls onto the BCG. The observed asymmetry is likely due 
to the cooling fiow being channeled along magnetic field lines. Based on the X-ray and Ha 
correlations, as well as the high effective conductivity in clusters with Ha filaments, it is 
likely that this cool gas is conducting heat from the surrounding ICM in order to remain 
ionized. Several questions remain unanswered in this scenario, such as the cooling process 
for the hot gas below ^ 0.5 keV, and the exact ionization mechanism for the warm gas. 
Upcoming papers in this series will carefully address these and other issues. 
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Properties of Individual Clusters 

In this section, we comment on the clusters in our sample for which the Ha emission is 
interesting in terms of its flux (or lack thereof) or morphology. We discuss possible formation 
mechanisms and implications on current theories as they apply to each case. 

Abell 0085 

A large ring can be seen in the X-ray image south of the cluster core, resembling the 
outer edge of a bubble. Unlike Abell 2052, however, this X-ray edge has no YLa counterpart. 
This cluster is one of three which has asymmetric Ha emission on small scales only. Based 
on the high Ha/NUV ratio, it is possible that the gas could be ionized by a recent burst of 
star formation in the BCG center. Previous works (e.g. Durret et al. 2005) have suggested 
that this cluster has undergone intese merging activity, based on the temperature/metallicity 
maps from XMM-Newton. This scenario is consistent with the relatively high temperature 
and entropy that we measure, compared to clusters with more fllamentary Ha emission. 

Abell 0133 

A single, thin Ha fllament extends northeast from the center of this cluster for ^ 25 
kpc. This fllament is only barely detected, but the detection signiflcance is increased slightly 
by the exact coincidence with an X-ray fllament. This cluster appears to be very similar 
to Abell 1795 in both morphology (long, thin strands in a single direction) and ionization 
mechanism (NUV/Ha ratio is consistent with ongoing star formation). Based on XMM- 
Newton observations, Fujita et al. (2004) came to the conclusion that this X-ray fllament 
is likely due to a buoyant radio bubble, since it was the only scenario put forward that did 
not directly oppose their results. However, we show that the metallicity in the fllament is 
consistent with coming from larger radii, not smaller radii, therefore it is inconsistent with 
the buoyant radio bubble scenario. 

Ahell 0478 

The X-ray morphology and Ha morphology match up very well in this cluster, as seen 
in Figure [4] and also the unsharp masked X-ray image in Sanderson et al. (2005). There 
appears to be two nuclei in both the Ha and X-ray emission, while the continuum image 



^ 50 ^ 



is that of a normal, undisturbed elliptical galaxy. The available UV data for this cluster 
is quite poor, but it appears that there is a significant amount of on- going star formation 
in the center of this cluster. However, the combination of shallow data and high extinction 
(E(B-V)=0.52) makes it difficult to trust these data. Without a deeper UV image of the 
cluster core, it is difficult to determine the heating source. This cluster also has the highest 
Ha luminosity of any other in our sample and a significant amount of structure extending 
beyond the optical radius of the galaxy, making it one of the most intriguing clusters in our 
sample. 

Abell 04 96 

This cluster has arguably the most interesting Ha morphology of any cluster in our 
sample. There are at least 5 distinct filaments, with various shapes and directions. The two 
longest filaments run parallel to each other for ^ 12 kpc. Using unsharp mask HST images. 
Hatch et al. (2007) found significant structure in the dust which may be coincident with 
some of our detected YLa filaments. Dupke et al. (2007) also discuss the presence of a cold 
front, as can be clearly seen in Figure [4j which may be due to an off-center encounter with 
a massive dark matter halo roughly 0.5 Gyr ago. Regardless of the method of producing 
this morphology in the X-ray, we note that the Ha filaments do not extend beyond the cold 
front in this cluster (^ 15 kpc), despite the fact that the cooling radius is nearly 50 kpc. 

Abell 0780 (Hydra A) 

Wise et al. (2007) present a summary of the X-ray and radio properties of this cluster, 
showing the excellent correlation between the radio jets and the X-ray cavities. The arcing 
Ha filament that we detect north of the BCG appears to be spatially correlated with the 
radio jet. The NUV emission for this cluster is very bright, but it is unclear if it is has 
morphological similarities to the Ha emission without higher resolution UV imaging. A very 
bright star to the west of the BCG made the data reduction slightly more complicated for 
this cluster, leading to a lower S/N in the region with the southern filaments. 

Abell 1644 

An arcing trail of cool gas extends south from the core of this cluster. This morphology 
is likely the result of an off- axis encounter with a subcluster (Reiprich et al. 2004). We detect 
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a curved Ha filament coincident with the cool x-ray gas, as well as two perpendicular Ha 
strands to the west and north of the cluster core. These previously-undetected filaments 
have large length-width ratios, suggesting that turblence must be low or the magnetic field 
strength must be high in this cluster core. Additionally, we find very little evidence for an 
AGN in this cluster, in terms of the radio power. X-ray morphology and hard X-ray flux. 



Abell 1795 

A pair of thin lOOpc), intertwined filaments extend south from the cluster core for 
^ 50 kpc (MV09). These filaments appear to be heated by strands of young stars which 
lie along the full extent of both filaments. However, the elevated [Nil] /Ha ratios found in 
MV09 and the literature suggest that either the star formation is dominated by massive 
stars, or another ionizing mechanism is also important. These twin filaments are coincident 
with a giant X-ray filament which is unresolved by Chandra. Additionally, this cluster has 
powerful radio jets which are offset in angle from the filament direction. The cooling wake 
hypothesis may be valid in this case, as the BCG appears to have a strong peculiar velocity 
relative to the cluster (Oegerle & Hill 2001). However, we note that a cooling wake should 
intuitively have slightly more turbulence than a straight cooling flow, which is inconsistent 
with the picture of long, thin filaments that we observe. 



Ahell 2029 

This cluster is the most massive in our sample and also has the most massive BCG. 
However, despite being exceptionally bright at X-ray, NUV and 1.4 GHz, we detect no Ha 
emission. A clue as to why this may be lies in the unsharp mask image of the X-ray gas. 
Unlike most other massive clusters, Abell 2029 is quite symmetric with no signs of any 
structure in the X-ray image. This suggests that the X-ray structure which we observe in 
all Ha-bright clusters is necessary to produce the filamentary Ha emission. This provides 
yet another link between the Ha filaments and the cooling X-ray gas and, at the same time, 
offers a counter-argument to the idea of buoyant radio bubbles being responsible for the 
emission. 
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Abell 2052 

The Ha emission surrounding the BCG in this cluster is coincident with radio-blown 
bubbles in the central region of the cluster. These bubbles to the north and south of the 
cluster core are filled with radio emission, which likely originated from the AGN within the 
BCG (Blanton et al. 2003). Since the Ha emission is seen primarily along the edges of the 
northern bubble, we suspect that shocks may be responsible for the heating in this case. The 
very low NUV/Ha ratio seen in Figure |6] further supports this hypothesis. As this cluster is 
more strongly infiuenced by the AGN than any other cluster in our sample, it is an outlier in 



some of our plots. As an example, in Figure 12 Abell 2052 is represented by the point with 



the lowest on-filament metal abundance. This abundance measurement is most consistent 
with the core abundance, not the abundance at the cooling radius. Thus it is possible that 
our simple picture of the Ha filaments being tied to the cooling fiow may not apply to this 
cluster. 



Ahell 4059 

The core of this cluster contains two peaks, only one of which is reproduced in Ha. 
The Ha emission extends for a short distance to the west, coincident with an extended arm 
in the X-ray image (Choi et al. 2004). The radio emission from this cluster extends north, 
perpendicular to the direction of Ha emission, suggesting that the radio jets are not playing 
a role in the heating. 



Sersic 159-03 

The optical emission in this cluster has been observed by Jaffe et al. (2005) and more 
recently by Oonk et al. (2010). Both authors find molecular gas which traces the ionized 
gas out to 20 kpc. We are able to detect Ha to 35 kpc, where the surface brightness blends 
into the background. It is likely that a deeper exposure would detect the ionized gas out 
to the cooling radius of ^ 60 kpc. There is no evidence for an AGN in the core of this 
cluster, despite the fact that it has some of the most extended and complex filaments. This 
suggests that the phenomenon of warm ionized filaments in clusters is likely independent of 
the presence of an AGN. 



